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Abstract 


Measurement of noise is important m engineering A method to measure power 
spectral density (PSD) of low frequency noise ( 1 / | / ) has been developed 

In this method “Constant Q” filtering of noise is done For ‘ constant Q” filtering 
purpose Wavelet transform is used Wavelet transform has perfect reconstruction 
and “Constant Q” properties All cases of noises (different 7 s ) are studied and 
algorithm is developed to give power spectral density of any type of noise Perfect 
reconstruction property of Wavelet transform is used to synthesize 1/ ] / |^ noise 
from white noise 
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Chapter 1 


Introduction 


Study of noise is of importance in science and engineering because it sets lower limits 
to the accuracy of any measurement In order to evalute these limits one must know 
the magnitude and sources of noise involved and must learn how to measure and to 
locate these Generally one uses two prototypes of noise (1) White noise and (2) 
flicker noise or coloured noise 

An easy way to distinguish between these noises is to consider their autocorre 
lation function Rxx{t) as given by Keshner [1] and shown in fig 1 1 



The white noise has the well known delta distributed autocorrelation function If 
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we represent the power spectral density of the noise as proprortional to 1 / | / I'’' then 
the white noise case corresponds to 7 = 0 The case of 7 = 2 is also well known 
because a variety of processes give rise to exponentially decaying autocorrelation 
function, examples of these processes are random telegraph signals, generation 
recombination noise in semiconductor etc The case of 7 = 1 is called ‘ 1/f noise ’ or 
‘flicker noise’ and has the largest autocorrelation as evident from fig 1 1 Moreover 
in engineering practice, one finds that ‘ 1/f noise is ubiquitous and appears m 
several physical phenomena Some examples of these are [1],[2] 

1 noise in the currents and voltages of electronic devices (eg field effect and 
bipolar transistors, schottky zener and tunnel diodes ) 

2 resistance fluctuations in metal film, semiconductor films and contacts ger 
manium filaments in carbon and aqueous solution thermocells and concen 
trations cells 

3 burst errors m communication channels 

4 physiological time series such as instantaneous heart rate records for healthy 
patients EEG vaiiations under pleasing stimuli, and insulin uptake rate for 
diabetics 

5 biological time series such as voltage across nerve and synthetic membranes 

The study of flicker noise has been widely persued over almost six decades and 
even now several unresolved issues regarding it’s modelling and it s affect on mea 
surcment and communication are remaining The emphasis in the present work has 
been on the measurement issues and, particularly, an attempt has been made to 
exploit the recently developed wavelet theory [3] [8] for analysing and synthesizing 

l/f noise as developed by Wornell [2] Our contribution m this work is m designing 
the filters and verifying the algorithms which relate to a theorem on 1/f noise as 
proved by Wornell [2] But before presenting this work, let us consider some simple 
schemes for the study of noise 
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Measurement of white noise can be done by constant bandwith filters as shown 
in figure 1 2 (a) Here filters are of constant bandwidth, therefore power output of 
each filter will be the same But if noise is 1 / ( / ]'’' and 7 = 1 , then it’s measurement 
can be done by “Constant Q”or Constant relative bandwidth filters as shown in fig 
12 (b) 



al 


a2 


a3 


a4 


aS 



pi 


p2 


p3 


p4 


p5 


LCD CONSTANT BANDWIDTH HI TERS 
al=a2=a3=a4=a5 

a and p are power oiilpul 


( b) CONSTANT Q FILTERS 
pl=p2=p3=p‘l=p5 


Figure 1 2 Constant Bandwidth And Constant Q Filters 

Constant Q filters can be designed [14] through a number of existing methods 
other than that of wavelet based designing We can get Constant Q filter bank by 
designing each filter individually But this method will require a lot of computation 
Also with this method reconstruction will not be possible Whereas in Constant 
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Q filter designing thiough wavelets, only one sequence (which, should fulfil certain 
conditions) is requiied Also this method ensures perfect reconstruction 

In the present work “Constant Q” filters have been designed through wavelet 
theory These “Constant Q” filters are then used to give power spectral density of 
different noise sequences 

To generate 1// noise time senes from white noise ( in order to analyse it by 
Constant Q filters)we have used two methods 

• Through digital filters designed using wavlet theory 

• Through digital filter obtained from analog filter using matched Z transforma 
tion for analog to digital conversion Analog filter has a slope of — lOdb per 
decade 

In first method a number of white noise sequences, having same power but differ 
ent lengths are required, whereas in second case only one white noise sequence is 
required 

Because in practice we can’t have ideil white noise sequence, we have used Knuth s 
[11] white noise sequence generator whose sequential correlation (shown in plot m 
ch 5) IS very near to that of white noise Because of this fact, noise after filtering 
can not be 1 // m the whole range 0 to tt 

1 1 organization 

This thesis is divided into six chapters 

• In chapter 2 wavelet theory is discussed In this chapter decomposition and 
reconstruction algorithm for any sequence is given Decomposition algorithm 
decomposes a signal in to different sequences in which each sequence corre 
spends to some fiequency band of input sequence This therefore corresponds 
to Constant Q filteiing of sequence 

• In chapter 3, a theorem [2] is given, by which any type of noise can be synthe 
sized from white noise In this chapter results are drawn from this theorem , 
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to analyze any type of sequence through Constant Q filter bank 

• In chapter 4 digital filter(s) design is given to generate 1// noise Here both 
filter designing cases are discussed as stated earlier 

• chapter 5 has results and relevant graphs 

• In chapter 6 we have drawn conclusion/comparison based on above results 
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Chapter 2 


Wavelet Theory 


Wavelet theory provides a unified framework for a number of techniques which had 
been developed independently for various signal processing applications 

In fact wavelet theory covers quite a large area It treats both the continuous 
and discrete-time cases It provides very general techniques that caji be applied to 
many tasks in signal processing, and therefore has numerous potential applications 
In particular, the wavelet transform is of interest for the analysis of non stationary 
signals, because it provides an alternative to the classical STFT (short time Fourier 
transform) or gabor transform The basic difference between the two cases is in 
the choice of wmdow’^used In STFT a simple analysis window is used, whereas in 
wavelet transform short windows at higher frequencies and long windows at lower 
fiequencies are used This is in the spirit of so called “Constant Q” or constant 
relative bandwidth fiequency analysis 

For some applications it is desirable to see the WT as a signal decomposition 
on to a set of basis functions In fact, basis functions called WAVELETS always 
underlie the wavelet analysis They are obtained from a single prototype wavelet by 
dilations and contractions as well as shifts The prototype wavelet can be thought 
of as a bandpass filter, and the constant Q property of the other bandpass filters 
(wavelets) follows because they are scaled versions of the prototype 

^By window we mean a signal having a window like structure which has negligible amplitude 
outside certain range 
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There are several types of wavelet transforms, and depending upon the applica 
tion, one may be preferred to the others For a continuous input signal, the time 
and scale parameters can be continuous, leading to a continuous Wavelet Trans 
form(CWT) They may as well be discrete leading to a wavelet series expansion 
Finally the wavelet transform can be defined for discrete time signals leading to a 
discrete WT In the latter case it uses Multiresoluiion signal processig technique In 
this thesis we have used Multiresolution signal processing technique as a tool to get 
information about PSD of l/( / P noise, so we will emphasize on this method 

2 1 A Review Of Multiresolution Analysis And 
Orthonormal Wavelet Bases 

The idea of multiresolution analysis [3] is to write functions f as & limit of 
succesive approximations, each of which is a smoothed version of /, with more and 
more concentrated smoothing functions 

More precisely, a multiiesolution analysis consists of 

(1) a family of closed embedded subspaces Vm C m & 

C ^2 C Vi C Fo C F-i C V-2 (2 1) 

such that (2) 

n I'm = w.uTT = m (2 2) 

mGZ m^Z 

and (3) 

/€ K. ^ /(2) G K*-!, (2 3) 

moreover, there is a ^ € Vo such that, for all m ^ Z , the (f>mn constitute an 
unconditional basis for Fn, that is, (4) 

Vm = linear 3pan{(l>mn,n e Z} (2 4) 

^Vector space of all measurable square integrable one dimensional functions 
^Set of integers 
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and theie exit 0 < A < B <cx. such that, for all Cnnez € 

(2 5) 

n n n 

Here — n) Let denote the orthogonal projection onto 

It IS clear from 2 1, 2 2, that = /, for all / e The condition 

2 Sensures that Vm correspond to different scales, while the translational invarience 
/ € Vm—*f{ — 2'"n) € Vm for all n 6 Z is a consequence of 2 4 and 2 5 <j!>{x) 
is known as scaling function or father wavelet 
Defining 

Cmnif) = <KnJ> (26) 

Then Pmf = E<^mn(/)<^mn (2 7) 

n 

Pmf is the approximated version of / in Vm vector space Now again if we ap 
proximate f at Vm-i then the information loss between two approximations can be 
obtained from projection theorem which is equal to the orthogonal projection onto 
orthogonal complement of Vm in Vm-i i e defining Wm as the orthogonal comple 
ment, in Vm-i, of Vm, 

V,n-l = K.0Vrm , (2 8) 

Equivalantely, 

Wm=-QmL\R) With Qm = Pn,-l-Pm (2 9) 

It follows immediately that all the Wm are scaled versions of Wq 

feWm^ /(2”‘ )eWo (2 10) 

and that the Wm are orthogonal spaces which sum to L^{R) 

L\R) =^Wm (2 11) 

me2S 

Wavelet theory ensures (eq 2 1 eq 2 5) that there exist V* € Wq such that ipmn 
contitutes unconditional bases for Wm 
This completes our analysis 
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2 2 SUMMARY Multi resolution analysis(MRA) 

To start our analysis we choose satisfying 2 1 2 5, then we get from 4>{x) 

through following procedure Since (f> £ Vq <f>(2x — n) E VLi, also sinceVo C V-i 
hence (f) € VLi Thus there will exist C„ such that, 

9i(a:) = X;<?n<^(2x-n) (2 12) 

n 

Then we define 

= E(-l)"<^n+i<^(2a: + n) (2 13) 

n 

i/) IS called mother wavelet and V’mn are dilated(or contracted) and scaled versions 
of mother wavelet 

2 3 Wavelet based decomposition and reconstruc- 
tion algorithm 

Let us assume that we have a sequence C® We want to decompose this signal m 
levels correspondoing to different frequency bands of input sequence 

To do this we will use multiresolution analysis, i e we choose ^ satisfying 2 1 
2 5 After using this decomposition algirithm, sequences obtained will correspond 
to different frequency bands of input sequence In other words this algorithm will 
correspond to passing signal through Constant Q filters 

2 3 1 decomposition algorithm 

From C° we form a function f as follows 

/=E„C”?i„n.or 

/(*) = Sn - n), 

since Vb is linear span of set {4>on,n € 2’} hence / € Vo 
Also from our construction Vq = Vi®W\, hence we can write, 

f = Pif + Qif (214) 
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where > bases for Vi, and ipi^ —i bases for Wi hence from 2 6 and 2 7 we 
can write, 

Plf = EC'ln(/)<^ln , C'i„(/) =< (f>^^ f> (2 15) 

n 

Where i\/ is approximation of / in Vi vector space 
Using Slimier technique we can write for Qif 

Qif = =< / > (2 16) 

n 

Where (^i/ is difference of information between f and it’s approximated ■version Pif 
Writing cl for Cin(f) and for di„(/) 

From 2 15 

cl = < > 

cl = j2-'^^^(l>{2-^x~n)f{x)dx 

cl = 2~^^‘^4>{2~^x — n)Y^c°<l)okdx from above 

d k 

cl — ^cl f 2~^^^(f>{2~^x — n)(i){x — k)dx 
k d 

cl = f 4>{x/2)4{x — [k — 2n))dx 

k d 

ci = E/^(^-2n)c2 (2 17) 

k 

Where h{n) = 2~^l^ f <f){xf2)f{x — n)dx , is known as lowpass analysis filter Looking 
at eq 2 17 it is clear that cl from c” can be obtained by filtering it through h[—n) 
and then subsampling it by 2 Similarly for dl, from eq 2 16, 

dl = < fin f > 

dl = < fin fok > from above 

k 

dl = X! 4 < fin fok > 

k 

4 = Y^cl2-^P f f{x/2-n)(l>{x - k)dx 

k 

dl = Y^cl2~'^P ( f{xf2)<f>{x -{k - 2n))dx 
k 

4 = Y^clg{k-2n) (2 18) 

k 
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(a)DECOMPOSITION ALGORITHM 



(b)RECONSTRUCTION ALGORITHM 


Figure 2 1 Decomposition and reconstruction algorithm 

wheie g{n) = 2“^/^ f r/j(xf2)4(x — n)dx, is known as highpass analysis filter Looking 
at eq 2 18 and fig 2 1(a) it is clear that from c° can be obtained by filtering it 
through g{—n) and then subsampling it by 2 After getting c\ from c° we can again 
apply the same analysis on to get c* and i e we can get lower resolution signals 
from higher lesolution signals by repeating the analysis algorithm just described 
This completes the analysis part 
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2 3 2 Reconstruction algorithm 

In reconstruction part (see figure 2 1) we have details(cf^ ,je [1, j]) and one approx 
imated signal(c^) and we have to construct original signal c° 

We start with lowest resolution i e from resolution, and construct a function 
a as shown below, 

a = + (2 19) 

n n 

a 6 

Clearly a € Vj and b E Wj, hence a 6 Vj_i , and we can write a = Pj-if 

a = Pj-i/ = E from tq 2 1 (2 20) 


P,-if = E<K + T,<<l’o. (2 21 ) 

n n 


Also 



= 

— 


(2 22) 

c’-^ = 

^j—ik ^ fvom 2 20 


= 

< <^j-in, S > from 2 20 

k k 


= 

<f>]k > < <f>j-ln,i’3k > 

k k 


= 

j <f>{xl2)<f){x - (n - 2k))dx + 42 _i /2 J if[x/2)<i>{x 

— (n — 2k))dx 

c’-^ - 

‘'n 

Yj oih(n — 2k) + Y ^k9i.‘'^ ~ 2A:) 

k k 

(2 23) 


1 e we have got from and 

Thus recursively using above algorithm we can get Reconstruction algorithm 
IS shown m figure 2 1(b) in a schematic fashion 

This completes decomposition and reconstruction algorithm which decomposes 
signal in to different frequency band signals and given these different frequency band 
signals, oiiginal signal can be constructed 

To be more familiar with decomposition, reconstruction algorithms and detaul 
and approximated signals we look at fig 2 1, and consider frequency spectra of these 
signals 
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Consider decomposition algorithm in which from (which will have frequency 
spectrum from 0 to ir ), first we get and dj^, where is obtained by filtering c° 
through h( n) and then subsampling it by 2 is obtained by filtering c° through 
g( n) and then subsampling it by 2 Here h( n) is half band low pass filter(pass band 
IS 0 to 7r/2), and g( n) is half band high pass filter(pass band is ir/2 to tt) Hence 
cl will have information from 0 to 7r/2 of c“ and d^ will have information from 7r/2 
to r of c° Again and d^ are obtained from through similar procedure c\ will 
have information from 0 to 7r/2 of c^ and dl will have information from 7r/2 to tt 
of cl In other words will have information from 0 to 7r/4 of c° and d^ will have 
information from 7r/4 to 7r/2 of c° 

Hence after applying decomposition algorithm on a given signal (c°),/ times, we 
shall be having (I + 1) signals, m which / signals will be detail signals(from to 
dl) and one will be approximated signal(c|,) 

Here dl will correspond to x/2 to ir of c“ 

Here d^ will correspond to Tr/4 to 7r/2 of 
Here will correspond to ir/8 to ir/4 of c° 

Here d^ will correspond to 7r/2^“^ to ir/2‘ of c° 

And c(j will coiiespond to 0 to 7r/2‘ of c° 

Hence we have divided original signal m to (1+1) signals correspondig to different 
frequency bands of original signal Moreover what we have achieved through this 
decomposition is good fiequency resolution at low frequencies and poor frequency 
resolution at higher frequencies (That is why it is called multi resolution analysis) 
We can see two graphs given in chapter 5 right now in which White and colour 
noises are decomposed m to different d„’s and c„ by the filters designed m chapter 
4 Looking at c* and d* (mentioned m these graphes as 1/2 resolution signals) one 
can Very easi ly see that in case of white noise both cl and d^ are rich m magnitude, 
showing that c® has good spectium over the entire range 0 to tt Whereas in colour 
noise case d,^ is very poor in magnitude as compared to cl, showing that c° has very 
small spectrum in ;r/2 to ;r, which indeed is the fact 
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We have discussed types of WT in the beginmg of this chapter There was a 
wavelet series expansion represantation of signal We will discuss this representation 
in next section 


2 4 Wavelet series expansion representation 

An orthonormal wavelet transformation of a signal /(s) is generally described in 
terms of synthesis/analysis equations [2] 

= (2 24) 

m n 

where 

C = <V'mn,/> (2 25) 

J f 

Considering synthesis equation i e eq 2 25 we can write 

<C = / /(a:)2-'"/V(2-'"x - n)dx (2 26) 

putting 2~'^x — n = 2~'^y, 

<C = //(si + 2”-r.)2-™/V(2-”!,)dy 
= //(si + 2'"n)((>„o(si)<isi 

<C = {/(i')»s/>mo(-ai)}|„2™„ (2 27) 

This representation is known as octave band filter band interpretation Thus for each 
m, dn can be obtained from filter and sample operation This representation is shown 
m figure 2 2 

Considering analysis equation i e eq 2 24, it can be viewed as multirate modu 
lation as depicted in figure 2 2 

To undeistand how we get multirate modulation diagram as shown in fig 2 2 
consider d^, from eq 2 27, it is clear that for each m, d^ correspond to samples with 
consecutive samples 2”* seconds apart 
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(a)OCrAVE BAND 
FILTER INTERPRETATION 


(b)MULTIRATE MODULATION 


Figure 2 2 Octave band filter bank and multiiate modulation interpretation 
Hence will have form as shown in fig 2 3 And can be written as 

<C=EO("-’‘o) P 28 ) 

no 

where n = 2~'^x from octave band filter representation, theiefore 

no 

no 

no 

15 















Figure 2 3 Sequence obtained from filter and sample operation 
C*V-mo(x) = ^CV>m0(®-2”‘n) 

n 

= EC2'’"^V(2'’"(a:-2"'n)) 

n 

= 5];d;;‘2-’”/V(2~’”x-n) 

n 

* V’mo(®) — 53^n^»nn(®) 

n 

or we can write 

E(C*«x)) = (2 29 

m m n 

= /(^) (2 30 

1 e synthesis equation can be viewed as multirate modulation as shown m figure 2 2 



Chapter 3 


Synthesis And Analysis of Low 
Frequency Noise 


In chapter 2 we have discussed decomposition algorithm, which decomposes a signal 
in to different frequency bands This algorithm corresponds to Constant Q filtering 
of input signal 

Now we will state a theorem [2] by which any type of low frequency noise can be 
synthesized through reconstiuction algorithm of chapter 2 This theorem also helps 
in making the choice of wavelet basis that should be used to synthesize a particuler 
noise Initially we will motivate the structure of this theorem and then a rigereous 
statement will be made 

In chapter 2 we discussed wavelet based series expansion Let us write the 
function f{x) as follows 

= ( 31 ) 

m n 

Now in order to represent a process, having power spectral density of the form 
(T^/ I / using equation 3 1 constraints have to be imposed on and ^^{x) These 
constraints aie 

1 should be a collection of mutually uncorrelated zeio mean random variables 
with variances [2] 


17 


i;ar{d”} (3 2) 

Here is length of sequence Note that = 2"”^ from fig 2 3 therefore 
var{cf™} = 0-22'"'!' 

2 Orthonorrral wavbelet basis i/)(a:) should be having at least 7/2 order regularity^ 
Now we will state the theorem [2] 

THEOREM 3 1 Consider any orthonormal wavelet basis with order regu- 
larity for some R> 1 Then the random process constructed via the expansion, 
f{^) - Em En CV'mnCx) 

where are a collection of mutullay uncorrelated zero mean random variables 
with variances, 

var {d^} = , 

for some parameter 2 R > 7 > 0 has a time averaged spectrum 
5,H=cr’E„2’"|V(2”u) I" 
that IS nearly Iffie 

< Sf{oj) < |5|V 

for some 0 < ctl^ < <oc , and has octave spaced ripple i e for any integer 
A 

I w |’^5'/(w) = I 2-^w |'’'5/(2'^u;) 

With the help of this theorem, an algorithm is developed m the present work to 
synthesize and analyze noises for different /s 

Before proceeding towards next section we simplify eq 3 2 to get an equation 
for power(P„) of sequence as follows 

var{d^]=cT^{llL^]\ (3 3) 

Note that 

var{d:;] = ^ZKV 

„=1 

^see Appendix A 
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Hence power(P„i), 


L/m 

= = (3 5 ) 

n=l 

or we can write, 

Pm = (3 6) 

3 1 Synthesis of 1/ | / noise 

To synthesize 1/ | / ['' noise for given 7 , frequency range and output sequence 
length Lou(, steps that should be followed are given in fig 3 1 Frequency range 
should be of the form tt to 7 r/ 2 *, where k & Z Also Lout should be a multiple of 2 
To understand synthesis more clearly we consider an example Suppose we want 
a sequence of length 1024 (It is a multiple of 2) with PSD tr^/ [ / p over frequency 
range tt to jt/S In other words we have specified the values of a^, 7 , frequency 
range and output sequence length, and using above theorem we want to construct 
a sequence, (from white noise sequences) having PSD very near to as desired 
For this example fc = 3 (see fig 3 1) and hence 

1 Three detail signals are required 

(a) 7r/2 to TT signal 

(b) 7r/4 to 7r/2 signal 

(c) tt/S to 7r/4 signal 

2 Length of tt/S to 7r/4 signal is 
Loui/2'' = 1024/8 = 128 

Since Lm = 2""* therefore m = -7 for ir/S to 7 r /4 signal Hence this signal 
will be represented as d~^ 

Similarly length of ir/4 to 7r/2 ««) signal = 256 
And length of t/8 to 7r/4 (d~®) signal = 512 

3 Power required of d'"” signal = P_r = from eq 3 6 
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* The number of detail signals 
IS obtained as {ntoTT/L, 

This fixes the number of bl 
ocks m fig 2 1(b) 








The leftmost block of fig 
2 1(b) corresponds to freq 


uency range TT ia ILand the 

iKl 

rightmost block corresponds 


to frequency rang e ^ 77^ 

The length of detail signal 
for ^/ 3 ^^ IS Louty^K 

note that ^°^where m 
represents scale factor, and 

denotes desired seq 
uence length of output 


This step calculates the va 
nance of the white noise 
asociated with fig 2 1(b) 
using eq 7 ^ and Leidr 

The sequence length of 
IS the same as of e:^,^ut the 

power of C^s zero Now all 

the inputs for fig 2 1(b) 
aie specified The upscaling 
opeiation and filter design 
aie given m chapter 4 


Figure 3 1 Steps foi aiianging the power and sequence length of input and c™ 
in fig 2 1(b) 


20 






Power lequired of cZ„® signal = P _8 = 

Power required of d~^ signal = P.g = 

Hence by using df“^, d~® and we can get through reconstruction algorithm 
which will have PSD very near to a^j 1 / m the range t /8 to tt 

Fronn above we have concluded following 

1 To synthesize 1/ \ f \'* noise for 1 > 7 > 0, lower resolution signals (higher val 
ues of m) must have lesser power than higher resolution signals (see eq 3 6 ) 
In other words to synthesize 1/ | / ]'' noise for 7 < 1 case we require zero 
mean uncorrelated sequences of different lengths As length of sequence in 
creases(resolution increases) power of sequence must increase 

2 To synthesize 1/ | / 1'^' noise for 7 = 1 case, for all resolution signals, power 
required is same In other words to synthesize 1/ | / ^ noise we require zero 
mean uncorrelated sequences of different lengths but of same power 

3 To synthesize 1/ | / l’^ noise for 1 < 7 , lower resolution signals (higher values 
of m) must have higher power than higher resolution signals (see eq 3 6 ) 
In other words to synthesize 1/ | / I'*’ noise for 7 > 1 case, we require zero 
mean uncorrelated sequences of different lengths As length of sequence in 
creases(resolution increases) power of sequence must decrease 


3 2 Analysis of 1/ | / |'>^ noise 

In analysis pait we are given a sequence and < 7 ^ and 7 of this sequence has to be 
found out This is done through “Constant Q” filter bank designed using wavelet 
transform (See fig 2 1(a)) Once we are given a sequence, and it has been passed 
through decomposition filters (“Constant Q” filters), we shall be having one approx 
imated and other detail signals These signals will be of varying length with each 
decrease in resolution (increase in value of m), length will decrease to approximately 
half of its pievious length Hence after passing given sequence through decomposi 
tion filters, we shall have Pm and Lm of each detail signal, and we want to know 
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and 7 From eq 3 5 we have 

Pm = a^l/LmV-^ (3 7 ) 

or 

Pm = + {I- 'y)\og2Lm (3 8) 

Looking at above equation, we can make following conclusions regarding 7 of input 
sequence from the knowledge of Pm , Lm and the fact that as resolution decreases 
length of detail signal decreases These conclusions d,re 

1 If power obtained in coarser detail signal is less than power obtained in finer 
detail signal, it is a 7 < 1 case 

2 If power obtained in each detail signal is same then its a 7 = 1 case 

3 If power obtained m coarser detail signal is more than power obtained in finer 
detail signal, it is a 7 > 1 case 

(*!*= We could reach to the same conclusions as stated above in 1,2, and 3 b} looking 
at fig 3 2 In this figure power spectra of 7 < 1 7=1 and 7 > 1 noises are shown 
It IS evident from this figure that power spectrum of 7 > 1 noise converges faster 
than 7 = 1 noise, whereas power spectrum of 7 < 1 noise converges slower than 
7 = 1 noise 

We are using “Constant Q” filters to analyse noise and if input noise has 7 = 1 
then power output of each of these filter will be the same But because of faster 
convergence of 7 > 1 noise at higher frequencies (higher resolution) power output 
of higher resolution signal will we lessei than power output of lower resolution sig 
nals Similarly because of slower convergence of 7 < 1 noise at higher frequencies 
power output of higher resolution signal will be higher than that of lower resolution 

signal **) 

Therefore by plotting a graph between log Lm and logp^ for different detail 
signals obtained from “Constant Q filter bank, we can find < 7 ^ and 7 In next 
section we have considered an example for analysis of 1 / ] / 1 noise 
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Figure 3 2 Powei spectia of 1/ ] / noise 

3 2 1 Example 

Let us assume that we have a sequence of length 1024, and we want to know 
value of and 7 over certain frequency lange Here frequency range will decide 
how many times decomposition algorithm has to be applied on the input signal If 
we want to know and 7 over fiequency range (tt to 7r/2''), then decomposition 

algorithm (shown in fig 2 1 (a)) will have h steps and there w ” be fc detail signals 
Let us assume that k = 4, hence first detail signal will be of length S 512 i e 
it can be represented as simileily next coarser detail signal can be represented 


4 - 
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as and so on Now we plot a graph between logj Z^and logj Pm, let the graph 
obtained be as shown m fig 3 3 



Figure 3 3 Plot between power and length for 7 > 1 case 

1 e coarser detail signals are having more power than finer detail signals, it then 
IS 7 > 1 case and value of 7 can be obtained from the slope of above curve, let slope 
be k, then 

7 — 1 = and value of 7 is obtained 

If graph obtained is as shown in figure 3 4 , then it’s a 7 = 1 case 
However if graph obtained is as shown in figure 3 5 , it is a 7 < 1 case and value 
of 7 IS 7 = 1 — fc, where k is slope of the plot 

After obtaining value of 7 from the slope, cr^ will be obtained from eq 3 8 
Hence from Constant Q filter bank ( realized by wavelet based decomposition 
algorithm ) we can get P S D of any type of noise 
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Chapter 4 


Digital filter design 


In [9], [10] Iff like spectrum is obtained but that doesn’t give an output time series, 
instead it first generates sequences through a recursive algorithm, then takes auto- 
correlation of each sequence and desired 1// spectrum is obtauned by talcing FFT 
of each autocorrelation and then averaging these FFTs Since autocorrelation is not 
a linear operation, we can not get a sequence having 1 // like spectrum from this 
method 

In chapter 3 we have discussed that Iff noise can be synthesized from white 
noise sequences obeying certain power constraints We will discuss design aspects 
of these filters in next section 

We will also generate 1 // noise sequence through digital filters from A/D trans 
formation using Matched Z transform, for analog to digital transformation with 
analog filter having a slope of lOdb/decade Designing of this digital filter is dis 
cussed m section filter 2 


4 1 filter 1 

To realize decomposition (“Constant Q”) and reconstruction filters, we have chosen 
Daubachies [3] 16 pt h{n), reproduced in table 1 , as low pass analysis filter This 
filter has very high regularity and obeys regularity constraint of theorem 3 1 for all 
practical cases of 7 s (2 > 7 > 0 ) 
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n 

h(n) 

g(n) 

0 

054415842243 

000117476784 

1 

312871590914 

000675449406 

2 

675630736297 

000391740373 

3 

585354683654 

004870352993 

4 

015829105256 

008746094047 

5 

284015542926 

013981027917 

6 

000472484574 

044088253931 

7 

128747426620 

017369301002 

8 

017369301002 

128747426620 

9 

044088253931 

000472484574 

10 

013981027917 

284015542926 

11 

008746094047 

015829105256 

12 

004870352993 

585354683654 

13 

000391740373 

675630736297 

14 

000675449406 

312871590914 

15 

000117476784 

054415842243 


Table 4 1 16 pt h(n) and g(n) used m decomposition and reconstruction algorithm 

Any other sequence can be chosen as h{n) if it obeys conditions 1 to 3 of theorem 
given in appendix A and regularity condition of theorem 3 1 

High pass analysis filter g{n) is obtained from h{n) by the relation [3], 

^(n) = (-l)"/i(-n + l) (41) 

If h{n) is of length 16 i e n € [0, 15] , n 6 Z, then g{n) will also be of length 
16 with n e [-14,1] , n e Z However g{n) is shifted to right so that the shifted 
version g'{n) will be zero for n ^ [0, 15] , n £ Z 

After choosing h{n) and g{n) as explained above, “Constant Q” filter bank is 
then realized by writting a computer program for cascade configuration shown in 
fig 2 1(a) 
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To realize reconstruction filter bank of fig 2 1(b) h{-n) and g{-~n) are used, 
however both are again shifted to origin There is one difference between real algo- 
rithm written on computer for reconstruction filter bank and fig 2 1(b) This is as 
follows In computer program a backward shift of 15^ is made after each su mmi ng 
block of fig 2 1(b) This backward shift is required because of delay caused by 
h{—n) and g{—n) filteis 


4 2 filter 2 

l/f can also be obtained from white noise by passing it through a filter having slope 
of lOdb/decade slope This filter (analog) is realized by pole zero model with poles 
and zeros are as follows 

polel=3 hz zerol=8 5hz 

pole2=30 hz zerol=95hz 

pole3=300 hz 

slope of this filter will be approx lOdb/decade Digital filter is obtained by 
Matched Z transform for analog to digital transformation with sampling rate of 
1000 hz (In Matched Z transform a zero is added at Z = — 1 if number of zeros are 
less than number of poles) Difference equation of this filter is 

y{n) = x(n) — 499x(n — 1) — 977x(n — 2)4- 526x(n — 3) 

4-1 961y(n - 1) - 1 087y(n - 2) 4- 123y(n - 3) (4 2) 

A digital filter satisfying above difference equation is realized by writting a computer 
program 


^Because filter length is 16 If filter length is I then shift should be / — 1 
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Figuie 4 1 Digital filter through A/D conversion 
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Chapter 5 


Results 


In this chapter, the input and output waveforms and power spectral densities of 
the matched Z and wavelet based filters are given For the input and output of 
the matched Z filter the “detail” and “approximate” decomposition are also given 
Furthermore, the synthesis of 1// noise through wavelet based filter bank has been 
presented 

5 1 The input white noise sequence 

An input white noise sequence has been used for both the filter 1 and filter 2 This 
sequence has been generated using Knuths[ll] algorithm and its autocorrelation 
and power spectral density (PSD) are shown in figs 5 1 and 5 2 respectively 
One notes that these are not ideal autocorrelation and ideal PSD for the sequence 
to qualify as white noise but are very good from the point of view of practical 
realization 


5 2 The output 1// noise of matched Z filter 

1 The PSD of the output of the matched Z filter designed in chapter 4 to 
provide 1// noise at the output, is shown m fig 5 3 Note that there is some 
deviation of this PSD from that of ideal 1// noise and this can be attributed 
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to the nonideality of the input white noise as shown m fig 5 2 

2 The time sequence (realization) of the input white noise and the output 1// 
noise are shown in fig 5 4 A perusal of these plots exhibits the difference 
between the appearence of the white noise and Iff noise as would be seen 
say, on an oscilloscope or a recorder 

3 A noise sequence can be described in three different ways, 

(a) by providing the realization, (b) by providing autocorrelation and PSD 
or (c) by providing the “detail and “approximate decomposition through 
wavelet analysis We have carried out the ‘detail” and “approximate’ decom 
position of all the waveforms studied m this thesis 

The white noise detail” and “approximate” decomposition is shown m fig 5 5 
to exhibit the features of multiresolution analysis based on wavelets Four blocks 
(i e four levels of resolution) of fig 2 1(a) have been used and the ‘detail and 
‘approximate” output of each block is shown m fig 5 5 The 1/2 ‘detail’ signal 
corresponds to the output of (vr — 7r/2) filter, and the 1/2 “approximate signal 
corresponds to the output of (0 — 7r/2) filter, these are shown in fig 5 5(a) by red 
and blue respectively The 1/4 “detail” corresponds to (7r/2 — 7r/4) filter output 
and 1/4 “approximate’ to (0 — 7r/4) filter output, and so on We will return to 

the considerations of this figure after presenting the ‘detail’ and approximate 
decomposition of 1// noise for comparison 

The “detail and approximate’ decomposition of 1 // noise (obtained from the 
matched Z filter) are shown in fig 5 6 It is interesting to compare the 1/2 
“approximate” and 1/2 detail’ decomposition of the white and 1// noise as given 
in fig 5 5(a) and 5 6(a) respectively One notes that the 1/2 detail” of 1// noise is 
very much smaller m magnitude than the 1/2 ‘approximate signal One also notes 
that whereas the waveform ratio’ of the ‘detail” to the approximate’ signals of 
white noise are the same for all the resolutions, this is not the case for 1// noise For 
1// noise the ‘detail” is comparatively small at 1/2 resolution and the ‘waveform 
ratio” of the ‘detail to the “approximate” signal incrases with decreasing (1/4 
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1/8, — > 1/16) icholution The lowei resolution coiiesponds to lower fiequency range 
in P S D , hence fig 5 6 shows that the “detail” signal mcrenses at lower frequencies 

The output 1// n ise of the wavelet based filter can also be analyzed m a similar 
fashion, hence is not presented heie 

5 3 Syntaesis of 1// noise through wavelets 

In the pievious ection the detail and “approximate” analysis of 1// noise obtained 
as output ot the matched z filter has been guen Now we obtarn 1// noise as the 
output of a wavelet based filter bank (filter 1) shown rn fig 2 1(b) The details of 
this filter are given in chapter 4 

The rnput white noise to this filter is the same as described earlier and is shown 
in fig 5 7(a) for convenience and comparison The output realization of 1// noise 
through this filter is shown in fig 5 8 

To conclude, we emphasize that whereas the wavelet based filters have the ca 
pabihty of both analyzing and synthesizing a required noise waveform, the matched 
Z filter have only the capability of synthesizing a noise waveform However, the 
wavelet based filter for synthesis requires the input white noise sequences of differ 
ent lengths but of same power, whereas the matched Z filter requires only one white 
noise sequence 
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Chapter 6 


Conclusion 


In the piebcnt woik mcasuiement of low fiequency noise through “constant Q” fil 
teung IS done All crscs of low fiequency noise are studied 

Simulation of low frequency noise his been done by two methods (1) Matched Z 
transform technique for analog to digital conversion with analog filter having slope 
of 10 db/decade, and (2) Wavelet transform based technique It has been found 
that these techniques give 1/ | / |'^ noise when white noise is applied to the input 
It has been found that time for measurement of noise through “Constant Q” 
filter IS more than that icquued m FFT method in which first autocorrelation of 
sequence and tlu n it’s FFT is taken (o get P S D of sequence However with this 
method modelling of 1/f noise is possible This is the greatest advantage of this 
method 
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Appendix A 


We know that decomposition and reconstt notion is achieved thiough filters h{n) and 
g(n) wheie 

h{n) = <<f>iQ,(i>on> (Al) 

g{n) = <i}iQ,(l>on> (A 2) 

Thus to constiuct h{n) and g(n) we have to start with satisfying 2 1, 2 2, 2 3, 
2 4, 2 5, and then we get h{n) and g(n) from A 1 and A 2 respectively However we 
can get h(n) and g(n) directly without bothering about <j> This is due to following 
theorem [3] 

THEOREM A 1 Let there be a sequence such that, 

1 I h{n) I I n I' <oc for some e > 0 

suppose also that mo(e) = 2“^/^ En h{n)e^'^ 

can be written as mo(e) = [|(1 + C‘'')]^[En/(^)®^"^] 

where 

4 En I /(»^) I I ” r ^ ^ ® 

5 Sup.enlE„/(ny'“l<2''‘‘ 
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Define 


g{n) = 

(-irA(-n + l) 

(A3) 

$(£) = 

(2T)->/='n^o{2-^£) 

J=1 

(A 4) 

il){x) = 

- n) 

(A 5) 


n 


Where ^{e) ts fourier transform of (f){x) 

Then '^x — n) define a multiresolution analysis, the f>jit are the 

associated oithonoimal wavelet basis 

Hence we can choose h{n) duectly satisfying first three properties of above theorem, 
and can apply multiresolution analysis on any signal 

A 1 Regularity 

A function f(x) is order regular if its fouiier transform F(w) decays according 
to 

/■(u;) ~ 3i(| a; 1“'^), w -+oc 

Wheie the notation 3^( ) is to be understood to mean that 
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